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Abstract Searching for radio pulsars typically requires a bespoke software pipeline to 
efficiently make new discoveries. In this paper we describe the search process, provide a 
tool for installing pulsar software, and give an example of a pulsar search. 
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1 INTRODUCTION 

Searching for radio pulsars is a complex and time-consuming task, however we are continually moti- 
vated to make further discoveries by the important and often unique science that flows from the study 
of these objects. The first pulsar discoveries were made by detecting individual pulses on pen-chart 
recorders, however even with large telescopes only the brightest and closest pulsars can be detected 
through this method. 

There are two aspects to the pulsar signal that distinguish it from other astronomical signals and 
man-made radio frequency interference (RFI). Firstly, the constant spin period of the pulsar leads to 
a very regular train of pulses with periods ranging from 0.001 to 10 s. Secondly, free electrons in the 
interstellar medium cause the pulse to be dispersed, causing a group delay at frequency v, 



where n e is the electron density along the line of sight to the pulsar. We term the path integral the 
dispersion measure, DM, and give it in units of cm _3 pc. Although this adds complexity to the analysis 
of pulsar data, it also provides a clear distinction between astronomical and terrestrial signals. Therefore, 
if we are to discover new pulsars, we must develop a strategy for detecting dispersed and highly periodic 
radio pulses. 

In this paper we will discuss the general pulsar search strategy (Section 1), describe how to install 
the required pulsar search software using PSRSOFT (Section 2) and give an example of how to process 
an observation from the Parkes Radio Telescope (Section 3). 

2 THE PULSAR SEARCH STRATEGY 

The search for radio pulsars is the search for periodic, dispersed, radio signals on the sky. Therefore, 
this defines the parameter space over which we must search: position on the sky, DM, and pulse period. 
Additionally, in order to increase sensitivity to binary pulsars for which the pulse period varies due 
to the Doppler effect, we may need to include additional search parameters such as acceleration and 
higher order terms. Luckily, these parameters are generally independent, and so we can form them into 
a hierarchical structure where we loop over each parameter in turn. Figure[T]shows a schematic overview 
of the pulsar search process, with the following steps: 

- Extraction — this covers obtaining the relevant data file from the telescope data archive, and con- 
verting the data to a compatible data format. 




(1) 
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- De-dispersion — here we change our single file with multiple frequency channels to a large number 
of time-series each correcting for a different value of DM. 

- Periodicity search — each of the de-dispersed time-series must be searched for periodic signals, 
or for bright individual pulses. 

- Filter candidates — periodicities detected in each of the time-series are grouped together to pro- 
duce the final candidate listing. 

- Optimisation — finally, the original data can be 'folded' to optimise the period and DM, and 
generate the diagnostic plots used to distinguish good candidate pulsars from RFI. 
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Fig. 1 Schematic overview of the pulsar search process. Data products are shown on the left 
and software processes are shown on the right. 



2.1 Data Extraction 

A typical pulsar search data file will be a recording of a single position on the sky, with 64 to 1024 
frequency channels sampled every 32 to 256 /is. Occasionally, multiple positions will be recorded into 
the same file, e.g. when using a multi-beam receiver such as that used in the Parkes Multibeam Pulsar 
Survey 1 1 , however even in this case software is provided to convert the data a single observation per 
file. Unfortunately different instruments often record data in different formats, and converting the data 
to a usable format may require some inside knowledge of the instrument. Some effort is being made, 
however, to adopt the PSRFITS data format across many different instruments, and large databases such 
as the Parkes Pulsar Data Archive have converted all the data to PSRFITS ®. 

2.2 De-Dispersion 

As the pulsar signals are dispersed by the interstellar medium, the pulse arrival time is delayed by dif- 
ferent amounts at each frequency and the delay is proportional to v~ 2 . To counteract this, the telescope 
signals are recorded in many narrow frequency channels, and the effects of dispersion are reversed dig- 
itally. This process, known as de-dispersion, is simply a matter of applying a correcting time offset to 
each frequency channel and summing the channels. However for search observations the amount of 
dispersion suffered by the potential source is unknown, and this process must be repeated with many 
trial values of DM. This produces a large number of single-channel timeseries, for which the subsequent 
processing can be performed independently. For large searches, it may be beneficial to distribute these 
time-series to many different compute cores, allowing for the search to be carried out in parallel. 
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2.3 Periodicity Search 

The Fourier transform is an ideal way to detect the constant periodic signals produced by the majority 
of known pulsars. The Fourier transform is a decomposition of a data series into a sum of complex 
exponentials (i.e. sine waves). In this way an infinite sinusoidal wave will form a single infinitely narrow 
spike at the frequency corresponding to the period of the signal. More generally any periodic signal will 
form a peak at the fundamental frequency, plus a series of harmonics, with magnitudes depending on 
the shape of the input signal. By transforming the input time-series in this way, any pulsar signals in the 
data should be among the strongest signals in the Fourier spectrum, even if they are undetectable in the 
original time-series. 

Whilst the fundamental frequency of the pulsar is the strongest signal produced in the Fourier spec- 
trum, pulsars with narrow pulses will have a lot of power 'lost' in the harmonics. By using a technique 
known as harmonic summing, it is possible to regain some of this 'lost' power. This is achieved by 
stretching the power spectrum by a factor of two, and summing it with the un-stretched power spec- 
trum. In this way, the power from the second harmonic will be added on top of the fundamental power, 
increasing the signal by more than the V2 increase in noise (The increase depends on the pulse shape). 
One can then select the best signals from this new spectrum. This can be repeated to add more har- 
monics. However since the noise is also increased upon every addition, the S/N will reach a maximum 
at a harmonic fold related to the pulse width. Therefore it is important to keep the results from each 
harmonic fold used, and the harmonic fold that a pulsar was detected strongest in can give an indication 
of the width of the pulsations. Since the harmonic summing uses the power spectrum, it is possible that 
the process will boost the S/N for signals that are not truly harmonically related. Alternative methods of 
coherently summing harmonics can be used to identify the signals that are tmely harmonically related. 

The principal algorithm used to perform Fourier domain searches is the Fast Fourier Transform El. 
This algorithm computes a discrete Fourier transform in the order of N log 2 N operations, for transform 
length N, The improvement over a traditional DFT algorithm is immense, especially when the number 
of data points is large. For example, with the number of data points at 2 23 , a typical observation length, 
the computation time is reduced by a factor of 105. Recent work has looked to accelerate this processing 
even further through the use of the powerful graphics processing units of modern computers. By using 
the highly parallel processing architecture of these chips, the time for all parts of the search may be 
reduced dramatically. 

2.3.1 Binary Motion Compensation 

Binary pulsars have an apparent period drift due to the Doppler effect. If the binary period is of the same 
order as the observation time, the changing period will mask the signal in the Fourier domain. This is 
due to the power being distributed in different frequency bins. 

These effects can be effectively removed from a time-series before processing, however to com- 
pletely remove all the effects the binary parameters must be known. There are unfortunately too many 
parameters to perform a search for all possible orbits in a reasonable time scale. By expanding the or- 
bital motion as a polynomial, it is possible to approximate the orbit by using acceleration and jerl£] By 
using this approximation, searching for binary pulsars can be done with only one or two extra parame- 
ters. There are three typical methods for removing this acceleration from the search process, either by 
re-sampling the time-series EE, using incoherent methods or by using matched filters in the Fourier 
domain ^. 

The 'simplest' method is to re-sample the time series with a new sample rate that varies with time, 
corresponding to a particular value of acceleration and jerk. This effectively reverses out the binary 
motion and so the corrected time-series can then be searched using exactly the same methods as the un- 
corrected time-series. Since the optimal acceleration and jerk are not known beforehand, it is necessary 
to search over a large range of parameters to get the optimal S/N. 



Jerk is used in this document as a convenient contraction for 'rate of change of acceleration' . 
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Alternatively, one can search for acceleration in the frequency domain. The number of bins which 
the pulse power is spread across for a given orbital acceleration is proportional to the observation time. 
Therefore this spread can be reduced by splitting the observation into n smaller segments. The signal 
strength in each part is however reduced by a factor of 1/n, however by adding the power-spectra along 
lines of constant acceleration we can recover l/y/n of the original signal strength ratio. Although this 
incoherent method is less sensitive, it is very fast. 

In order to recover all the original signal strength in the frequency domain, we can use matched 
filters. When using matched filters, the time-series is Fourier transformed as normal, however the result- 
ing spectra are then analysed by convolving with the known response from each trial acceleration and 
jerk. This methods does not require repeated Fourier transformations, and so can be quicker than the 
time domain methods. Additionally, since this method searches over a parameter space in the frequency 
domain, where the search is actually done, the parameter space is never over-sampled. 

2.3.2 Single Pulse Search 

The flux density of individual pulses can vary considerably, and it is possible to detect many pulsars 
by the detection of individual strong pulses. In addition to this, there are classes of pulsar-like objects 
that have been shown to show so few total pulses that standard periodicity searches cannot detect them 
®. The detection of single pulses relies on the pulses being highly dispersed relative to terrestrial burst 
emission. Pulsars can be detected by one or more dispersed pulses that are significant above the general 
noise, and not associated with any un-dispersed terrestrial interference^. 

2.4 Candidate Filtering 

The periodicity search typically detects many hundreds of signals for each trial DM. Fortunately, many 
of these are harmonics of stronger signals, and many signals are detected at multiple DM trials. To 
reduce this to a manageable number of candidates, the signals are grouped together by period, keeping 
information about which DM trials it was detected in. This method can typically reduce the number of 
candidates to ~ 100 per observation. 

2.5 Time Domain Optimisation 

Once we have a small number of candidates, we can return to the original timeseries and 'fold' the data 
at the candidate period. This process steps through the time-series adding each sample to an appropriate 
output bin depending on the computed pulse phase at that moment. To fold data correctly requires more 
accuracy in the period than is provided by the finite width bins of a Fourier analysis. Therefore it is 
common to split the time-series into a number of sub-integrations, and fold those separately. With the 
shorter folds, the pulse will drift less in phase and therefore have less broadening. This also allows for 
the measurement of the variation of signal over the observation, providing useful clues for identifying 
pulsars. 

Summing the sub-integrations is equivalent to folding the entire time-series. Shifting each sub- 
integration in phase before summing can compensate for errors in the initial folding period. By shifting 
and adding the sub-integrations in this way, one can perform an approximate folding of the time-series at 
a number of different periods in a computationally efficient manner. In a similar fashion, using quadratic 
or cubic shifts with time can compensate for orbital acceleration and jerk. When the offsets correspond 
to the correct pulse period, etc. the power from the pulse should all fall in a minimum set of bins. The 
correction of the sub-integrations increases the effective signal strength, and therefore the S/N can be 
maximised. The way that the signal strength varies with period shows how well aligned in phase the 
pulses are, and so this information is useful when distinguishing highly accurate pulsar signals from 
man-made interference. 
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2.6 Candidate Selection 

Finally, we have produced candidate pulsars ready for inspection. Unfortunately, pulsar search pipelines 
often produce hundreds of candidates per observation, which often precludes manual inspection of each 
candidate. A simple method of selecting candidates for evaluation is to sort all candidates by S/N, 
however this often returns many instances of man-made RFI. More recently, graphical selection tools 
ESI or artificial neural nets ED have been employed to allow for a more efficient selection of pulsar 
candidates. 

3 INSTALLING SOFTWARE WITH PSRSoft 

Although the pulsar search process is conceptually simple, it can take a lot of time and effort to cre- 
ate a efficient and effective pipeline. Fortunately, many tools have already been developed for pulsar 
searching, and new algorithms are frequently being developed. A typical search pipeline may, however, 
depend on many different software packages. To help with installing this large variety of software, a 
software package installation tool PSRSOFT is provided. 

This section gives an overview of how to set up a pulsar search system with PSRSOFT. Before you 
install PSRSOFT you must ensure that you have a PYTHON interpreter and a C, C++ and FORTRAN 
compiler. PSRSOFT is tested using the GNU compiler collection (GCQ^]), however other compilers may 
work. It is essential that your C and FORTRAN compilers are compatible and can cross-compile, this 
typically means that they must be the exact same version. Some packages also require an Oracle/Sun 
compatible Java Development Kit. 

3.1 Installing and configuring PSRSoft 

PSRSoft is available from |http : / /www .puis arastronomy . net /pulsar/ sof tware/psrsof t| 
and is packaged as a gzipecj^TAP^Jfile. Installing PSRSoft simply requires unpacking the TAR file, 
typically using a command such as: 

$ tar -xzf psrsof t . tar . gz 
$ cd psrsoft 
$ Is 
bin/ config/ 

3.1.1 Configuring PSRSOFT 

Let us assume that you have unpacked the PSRSOFT install file as above. PSRSOFT is configured 
through a settings file, conf ig/prof ile, which defines variables using the BASH syntax. An ex- 
ample file is included which contains default settings that will likely work for your system. Begin by 
copying the example file: 

$ cp . /conf ig/prof ile . example . /conf ig/prof ile 

Most users will not need to change any settings in the profile file, however you can open it with any 
text editor to modify settings such as the preferred C or Fortran compilers. It is important to note that 
the C and Fortran compilers must be of the same version, typically gcc and gf ortran. You can check 
that the software versions are the same by confirming that the output of the following two commands is 
identical. 

$ gcc -v 
Using built-in specs . 
Target: i4 8 6-linux-gnu 



2 See: http://gcc.gnu.org 

3 See: http : //www . gzip . org/ 

4 See: http : / / www . gnu . org/ sof tware/tar / 
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Configured with: . . /src/conf igure ... 
Thread model: posix 

gcc version 4.3.2 (Debian 4.3.2-1.1) 

$ gfortran -v 
Using built-in specs . 
Target: i4 8 6-linux-gnu 

Configured with: .. /src/conf igure ... 
Thread model: posix 

gcc version 4.3.2 (Debian 4.3.2-1.1) 
3.2 Building Pulsar Software with PSRSoft 

PSRSoft is invoked with the command . /bin/psrsof t, taking a list of program names to install as 
arguments. A complete listing of software that are available can be found using the — search option. 

$ . /bin/psrsof t — search 
==== PSRSOFT version 1.5 ==== 
Pkg Index: 'stable' 



Updating package index 

Getting latest package descriptions... done 
Searching for package . * in stable tree 

1) psrchive 13.4+ 

2) cfitsio 3090 

3) f ftw 3.1.2 

4) presto-core 1.01 



If one piece of software relies on another, PSRSOFT will automatically download and install the 
dependencies before installing the requested software. To install all the software used in the example in 
Section 3, we can run the following. 

. /bin/psrsoft sixproc dspsr psrchive pulsarhunter 
==== PSRSOFT version 1.5 ==== 
Pkg Index: 'stable' 



Analysing dependancies 
Packages to be installed. . . 



1 


cfitsio 


NO 


(3090 2009-10-06 14:47) 


2 


f ftw 


N 


(3.1.2 2011-01-18 14:39) 


3 


pgplot 


N 


(5.2 2012-02-10 09:37) 


4 


sixproc 


N 


(5.0.5 2011-04-07 09:49) 


5 


tempo2 


NO 


(2012.1 2012-02-10 09:45) 


6 


tempo 


N 


(11.010 sf 2011-02-12 12:30) 


7 


psrxml 


NO 


(1.05 2010-02-04 12:30) 


8 


psrchive 


N 


(13.4+ 2012-04-11 11:01) 


9 


dspsr 


N 


(2.0 2012-04-12 14:52) 


10 


pulsarhunter 


N 


(1.3r79 2011-02-23 15:21) 


Install 


10 packages 


into psrsoft/usr? (y/n) 



The 10 packages are made up of the four packages we requested, plus the dependencies. The char- 
acter N indicates that the package is new (U would indicate that the package needs to be updated) and 
the symbol implies that the package is optional. For instance, if you do not require PSRFITS support, 
you can use the — no-cf itsio option to skip building cfitsio. Otherwise, answering 'y' to the 
prompt will begin the download and install of the requested packages. 
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4 EXAMPLE PULSAR SEARCH 

For this example we will download the file S13405_l.sf, which has source name G123070 from the 
Parkes Data Archive 12 '. Although this is a normal survey observation, it was pointed close to the very 
bright pulsar PSR J0835-4510 which should be easily detected by our search. 

4.1 Extract Data 

The downloaded data normally need to be converted to a format that is understood by the search software 
that is to be used. We will use the f ilterbank command from SlGPROC0to covert the downloaded 
data into a SIGPROC compatible format. 

$ tar -xvf PulsarObs . tar 
atnf-attribution . txt 
S13405_l.sf 

$ filterbank S13405_l.sf > S13405_l.fil 

We can check the details of our data by using the header command. 

$ header S13405_l.fil 
Data file 

Header size (bytes) 
Data size (bytes) 
Data type 
Telescope 

Datataking Machine 
Source Name 
Source RA (J2000) 
Source DEC (J2000) 
Frequency of channel 1 (MHz) 
Channel bandwidth (MHz) 
Number of channels 
Number of beams 
Beam number 

Time stamp of first sample (MJD) 
Gregorian date (YYYY/MM/DD) 
Sample time (us) 
Number of samples 
Observation length (minutes) 
Number of bits per sample 
Number of IFs 

4.2 Dedispersion 

Dedispersion is performed using the dedisperse_all command. This dedisperses to a large number 
of DM trials at once, with the DM range specified with the -d option. Since the input data file may be 
large, the software reads a block of data at a time, the size of which is specified with the -g option. Here 
we dedisperse from to 500 cm~ 3 pc using a block size of 100000 samples. 

$ dedisperse_all S13405_l.fil -d 80 -g 100000 

This will produce a large number (~ 326) files with the ' .tim' extension, each of which is a single 
time-series de-dispersed at a particular value of DM. 



352 

17154048 

filterbank (topocentric ) 
Parkes 



G123070 

08:32:37.0 

-44:51:00.0 

452.062500 

-0 . 125000 

256 

1 

1 

48664 . 615312499998 

1992/02/12 

300.00000 

536064 

2 . 7 

1 

1 



5 Here we use a modified version of SIGPROC available from https://github.com/SixByNine/sigproc, or 
through the PSRSOFT package 'sixproc' 
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4.3 Periodicity Search 

We can continue to use SIGPROC to do the periodicity analysis with the seek command. This program 
reads in one of our time-series and performs the Fourier transform and harmonic summing, outputting a 
list of periodicities and S/Ns. We can run seek on a single tim file, 

$ seek S13405_l . f il . 0000 . 00 . tim -fftw -head 

which has output a .prd file that contains a short header followed by columns of S/N and period in 
milliseconds. Each pair of columns is a different harmonic fold, which by default is 1,2,4,8 and 16. The 
period of PSR J0835^510 is 89 ms, and the pulsar is so bright that even in the 0-DM timeseries it 
should be the strongest candidate with a S/N of greater than 30. 

In practice we can run seek on all our .tim files using a simple shell script. In BASH we can run 

$ for timfile in S13405_l . f il . * . tim ; do seek $timfile -fftw -head ; done 
$ cat S13405_l . f il . * .prd > S13405_l.prd 

which will process all of the files and concatenate the output into a single file that we can use to identify 
candidates to begin the optimisation process. 

4.4 Candidate Matching 

Our combined .prd file contains more than 10000 signals, which is clearly far too many to analyse 
further, however many have the same period at different DM trials, or are harmonics of other stronger 
signals. We can reduce the number of candidates to something more manageable by combining signals 
with the same period or harmonics of the same period. For this we can use the best command from 
SIGPROC, however PULSARHUNTER provides more advanced matching and output options through 
ph-best. 

$ ph-best S13405_l.prd S13405_l — maxresults 50 

$ head -n 3 S13405_l.lis 

S13405_l_001.phcx.gz 445.600000 89.29830545 69.324158 

S13405_l_002 .phcx.gz 61.700000 267.79276611 57.694286 

S13405_l_003.phcx.gz 44.500000 290.20468989 69.324158 

The main output of ph-best is a text file with the .lis extension. This has columns of: can- 
didate name; S/N; period (ms); DM (cm _3 pc); acceleration (ms -2 ); jerk (ms -3 ); number of harmon- 
ics matched; and best harmonic fold. The top candidate should have a S/N of around 450, period of 
~ 89.29 ms and DM of ~ 69 cm~ 3 pc. 

In addition, ph-best produces candidate files with the . phcx . gz extension, which contain all 
the relevant information about the candidate in an XML format. 

4.5 Folding and Optimisation 

Once we have reduced the number of candidates to a reasonable number, we can fold the signals to 
produce diagnostic plots useful for identifying the candidates that are suitable for follow-up at the tele- 
scope. We can fold the data using DSPSR on the original . f il file. DSPSR can fold many signals at 
once, however for this example we will fold each candidate in turn. We can specify the folding period, 
in seconds, using the -c option, and the DM using the -D option; we can obtain these values for each 
candidate from the .lis file. Other useful options are -L to specify a sub-integration time (in seconds), 
-t to fold using multiple threads, -e to specify the output extension and -U to limit the memory usage 
to a value small enough to fit within the CPU cache. 

$ dspsr S13405_l.fil -c 0.08929830545 -D 69.324158 -L 10 -t 4 -U 1 -e subint 
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This produces a number of .subint files, which we can add together using psradd from 
PSRCHIVE. It is also useful to reduce the number of frequency channels to, for instance, 8. 

$ psradd *. subint -o S13405_l_001 . ar 
$ rm * . subint 

$ pam — setnchn 8 -m S13405_l_001 . ar 

Now that we have a psrchive archive of our pulsar file, we can use any of the psrchive tools 
discussed in van Straten & Demorest (in prep). In particular we want to use pdmp to optimise the S/N of 
the candidate by fine-tuning the period and DM. This also produced the diagnostic plots that are useful 
for distinguishing pulsars from RFI and noise. The phcx . gz files we created can also be updated with 
the additional information from pdmp, however pdmp also produces a postscript output that can be 
quickly viewed. 

$ pdmp -input-phcx S134 05_1_00 1 . phcx . gz -output-phcx S134 05_1_0 01 . phcx S13405_l_001 . ar 



$ gzip -f S13405_l_001 .phcx 

Since we need to do this for each of our candidates, it is again useful to create a script to automate 
the process. Below is an example BASH script to fold and optimise each of our candidates. 

# ! /bin/bash 
f ile=S13405_l 

while read line ; do 
set — ${file} 

cand= 'basename $1 .phcx.gz 1 # strip file extension 
psriod =, echo $2 | awk "{print $1/1000.0}" # convert to seconds 
dm=$3 

dspsr ${file}.fil -c $ {period} -D ${dm} -L 10 -t 4 -U 1 -e subint 
psradd *. subint -o ${cand}.ar 
rm * . subint 

pam — setnchn 8 -m ${cand}.ar 

pdmp -input-phcx $ { cand} . phcx . gz -output-phcx ${cand}.phcx ${cand}.ar 
gzip -f ${cand}.phcx 
done < ${file}.lis 



4.6 Viewing Output 

The postscript output from pdmp can be viewed with any standard postscript viewer, however more 
information can be obtained from the .phcx.gz files. These files can be viewed using either the 
JReaper graphical candidate selection tool 1 or individually using the ph-view-phcx command. 

$ ph-view-phcx S13405_l_0 01 . phcx . gz 

or to make a png image, 

$ ph-view-phcx S134 05_1_0 01 . phcx . gz --imageoutput 
$ open S134 05_1_001 . phcx . gz . png 

The candidate plot for our example candidate S134 05_l_001 is shown in Figure [2] The key fea- 
tures are that there is a clear peak in the 'DM curve', i.e. S/N as a function of DM and the pulsar is visible 
in all (or most) sub-integrations and sub-bands. Identifying pulsars and distinguishing them from RFI 
can be a difficult task, however this becomes easier with experience. 



M. J. Keith: Installation and Use of Pulsar Search Software 



File: S13405 1,001. phcx.gz RA08:32:37.0 Dec:-44:5 1:00.0 GI:-96.99 Gb:-2.97 MJD:48664.61 
ObsFreg 426 1MHz Tobs'160.0s SourcelD:G123070 Telescope:UNKNOWN 
SpecSNR:445.6 ' 
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Fig. 2 Example candidate plot showing the bright pulsar PSR J0835-4510. Panels are, from 
top to bottom, on the left side: signal strength as a function of frequency channel and pulse 
phase; signal strength as a function of sub-integration and pulse phase; and signal strength as 
a function of pulse phase. On the right side: S/N as a function of trial period and DM; and 
S/N as a function of DM. 



If a good candidate pulsar is identified, it is recommended to ensure that it not already known by 
consulting the ATNF pulsar catalogu^] If the candidate is not a harmonic of a known pulsar, then con- 
firmation of the pulsar can be done with follow-up observations, or by searching previous observations 
of the same location in archival data. 

5 CONCLUSION 

Although there are many different ways to search for radio pulsars, most common techniques follow 
the same steps as described here. Utilising data from different telescopes can bring new challenges, 
however most data formats are relatively simple and can easily be converted. In future, it is likely that 
all telescopes will generate data in the same PSRFITS format as used by the Parkes Data Archive, which 
is readable by most pulsar search software. Therefore, by following the principals in this paper, you 
should be able to process any pulsar search data available. 
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Appendix A: EXAMPLE PROCESSING SCRIPT 

This appendix contains a BASH script that will process a provided PSRFITS format file. Readers can 
download from http://zmtt.bao.ac.cn/psr/soft/processPSRFITS.csh 
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